Interplay Between Structural Randomness, Composite Disorder, and Electrical 
Response: Resonances and Transient Delays in Complex Impedance Networks 
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We study the interplay between structural and conductivity (composite) disorder and the col- 
lective electrical response in random networks models. Translating the problem of time-dependent 
electrical response (resonance and transient relaxation) in binary random composite networks to the 
framework of generalized eigenvalues, we study and analyze the scaling behavior of the density of 
resonances in these structures. We found that by controlling the density of shortcuts (topological 
randomness) and/or the composite ratio of the binary links (conductivity disorder), one can effec- 
JH \ tively shape resonance landscapes, or suppress long transient delays in the corresponding random 

impedance networks. 

PACS numbers: 89.75.Hc, 84.30.Bv, 05.60.Cd 

' ^ ' • Resistor networks have been widely studied since the 70 's as models for conductivityproblems and classical transport 
(— I in disordered media [1] . With the recent surge of research on complex networks 0, Q , resistor-networks and related 
I flow models have been employed to study and explore community structures in social networks [3| and to construct 
• j recommendation models for community networks Also, resistor networks, as abstract models for network flows 
' ^ • with a fundamental conservation law Q, were utilized to study transport in scale- free (SF) 0, [1, P , in a class of 
self-similar in small- world (SW) networks and in tree structures and hierarchial lattices [ij 

Complex impedance networks have been investigated to study electrical and optical properties of two-dimensional 
thin films 1^ and dielectric resonances of two-dimensional regular lattice structures, lattice animals, and other 
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fractal clusters 15, 1^, 17, 18|. In this Letter, we investigate electrical response (resonances and delays) when both 
^ the structure and the composition of the local conductances can be random, and we focus on the interplay between 
^ 1 structural and composite disorder, and response. Random structures, in particular, random nanowire networks, can 
play a key role in the design and fabrication of future electronics devices, such as transistors or interconnects 3, 2^, 



2ll . |22| | . Assessing performance and reliability of these systems requires to understand the time-dependent intrinsic 
' ' ' electrical response (resonances and transient delays) of these devices which have to switch electric currents on and off, 
and driven by high clock speeds. For example, in a random nanowire network made of single-wall carbon nanotubes, 
the individual wires can be either conductors or semiconductors (based on their individual chiralities), resulting in links 
. with (binary) composite disorder; their natural composition comes with a dominance of the semiconducting tubes [2^ . 
' Likewise, inherent delays in electrical signal propagation can have crucial effects on processes in neuronal networks. 
■ The compartmental-model representation of passive dendritic trees is an R1C-R2 network (each compartment consists 
of an Ri membrane leakage resistor in parallel with a capacitance C, and compartments are connected with an R2 
0^ junctional resistor) [23]. The framework employed here can be employed to study the effects of local defects (damaged 
or destroyed links) on global signal delays, ultimately governed by the structure and link disorder in the network. 

Here, we focus on the interplay between topological randomness, conductivity disorder, and system response. While 
the resonance and relaxation properties are well understood in low-dimensional structures with conductivity (bond) 



disorder 17|, and recently on the complete graph |2J,[25|, to our knowledge, a similar investigation on complex random 



d network structures with link conductivity disorder have not been initiated or explored. We employ the framework 



applicable to binary link disorder [13|, [15|, [17|, |2J, |25|. The powerful feature of the framework is that it can be 
employed to study the singularities of the electrical response associated with any kind of binary link disorder on any 
graph (random L-C, RL-C, R-C, or more complicated composite circuits, involving two, but individually arbitrarily 
complicated building blocks). 

The equations governing current flows in any network can be written as (il. [ill [l8| '^jC^ijiVi — Vj) — I{5is — 6it), 
where aij are now the possibly complex link conductances (or admittances). Nodes s and t are the nodes where a 
current / enters and leaves the network, respectively. The above equation can be rewritten as '^jLijVj = I{Sis — Sn), 
where = Sij X^j^^i '^n ~ '^ij ^he Laplacian of the underlying graph with complex couplings, also referred to as 
the admittance matrix in the present context. 

For example, in L-C composite networks (or in RL-C composite networks with weak dissipation), in order to find 
resonance frequencies one can identify the non-trivial singularities, the "zeros" of the admittance matrix (corresponding 
to the "poles" of the complex impedance matrix), i.e., requiring that zero input current gives rise to finite potential 
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FIG. 1; Network structures: (a) one-dimensional ring; (b) small-world network with random shortcuts added between nodes 
on a ring (Id SW) [2^; (c) complete graph (each node is connected to all others). In binary composite networks, each existing 
hnk has a complex conductance ai or 02 with probability q and 1—q, respectively. 



differences in the network jlT], |l8| . One can also show that the zeros of the conductivity matrix are directly related 
to the transient relaxation times in R-C composite networks, characterizing how fast the system responds to step-like 
on/off signals. 15l. Il7|. For binary composite networks, the conductance disorder of existing links in the structure 
is characterized by a single parameter (composite ratio) g, such that Uij = ai with probability q, and aij = ai with 
probability (1 — q) (and obviously, Oij = if nodes i and j are not connected). For example, for an L-C composite 
network, a\ = iujC, (T2 = {iujL)~^, while for an R-C composite network, ai — iujC, (72 = R^- Hence, for resonance 
condition in L-C networks (and for relaxation times for R-C networks), one searches for the nontrivial solutions of 
'^jLij{uj)Vj — 0, or L(aj)V = in a more compact notation. Then for any fixed graph and any realization of 
the binary link disorder, one can rewrite the above expression for the resonance condition (or to extract transient 



relaxation times) 13, 2 
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Here, Hij = Sij J2i=£i ^ where /ly = -1, 



-1, if cFij— (7i, (72, 0, respectively. Similarly, F^ 



(1) 



is just the (topological) network Laplacian of the underlying graph (h~^ — 1,0 is obviously the adjacency matrix of 
the network) . The expression for A in Eq. ([T]) establishes the connection between the generalized eigenvalues and the 



resonance frequencies of LC, or the transient relaxation times r, of RC composite networks 
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Hence the above generalized eigenvalue problem Eq. ([T]), where T is real symmetric and nonnegative and H is real 
symmetric, provides a framework to identify the resonance frequencies (density of resonances in the large TV limit) 
or relaxation times in the respective binary composite networks. It is also clear from the above framework that the 
resonance (and relaxation) spectrum (except from pathological cases) is independent of the choice of nodes where the 



current enters and leaves the system, thus, they represent intrinsic characteristics of the network |15l . llTl 

In what follows, for brevity, we use the "resonance" terminology in composite networks, and will also refer to A as 
"frequency". The eigenvalues of the above system always fall in the [—1, +1] interval, Ai < A2 < . . . < Xn, with true 
resonances corresponding to — 1 < A^ < 1. We focused on two important observables, the density of resonances p{\) 



and number of resonances per node p 15,17, 13, 2^ 



P(A) 



p{X)dX 



riR 
N 



(3) 



where nn is the total number of true resonances (not associated with Xj = ±1). In this work, we determined 
the spectrum of the generalized eigenvalue problem Eq. ^ numerically, and constructed the above observables by 
averaging over 10,000 realizations (1,000 realizations for the largest system size) of both structural and composite 
disorder. 

Before studying random structures with binary composite disorder, we recall two known extreme cases: the one- 
dimensional ring [Fig.llfa)] and the complete graph [Fig. [Tl^c)], both with the same composite (binary link) disorder. 
For a one-dimensional ring, there is a single resonance frequency [l3| (which can also be obtained via elementary 
considerations). More specifically, nji ^ 1 and the frequency is distributed binomially about {X) = 2q — 1. Thus, 
in the large- A'^ limit, the density of resonances Eq. ^ approaches a Gaussian distribution with the above mean and 
vanishing width, i.e., a delta function. For example, for g=l/2, p{X) ~ ■^'^(A) and p : 
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FIG. 2: (a) Average number of resonances per node vs density of shortcuts p in a Id SW network with composite ratio q=l/2. 
The inset shows the same data on log-log scales, with the straight dashed line indicating the asymptotic large-A'' behavior, 
p ~ p. (b) Density of resonances for fixed number of nodes A^=1000 with q=l/2, and for different values of p. (c) Density of 
resonances in the large-p regime for 5=1/2, for various system sizes. The inset shows the scaled plot of the same data, p{X)/p^^^ 
vs p^^^X. These and the following plots all show ensemble averages over 10,000 network- and composite-disorder realizations 
(1,000 for the largest system size). 



In the other extreme case where all nodes are connected to all others (i.e., the complete graph), using a path- 



integral approach [2J, |25|, |26[ in the large-A^ limit, Fyodorov obtained (without loss of generality, for q=\/2) that 



p[X) ^JI^ ^(^) and p ^j:^ ^ J ^Yie total number of resonances approaches the number of nodes, but they are all 
narrowly centered about the same frequency (becoming fully degenerate as N—^oo). 

Now, we consider small-wor ld ( SW) networks [23] as random structures, where random shortcuts were added to a 
one-dimensional ring (Id SW) [Fig.IHb)], resulting in an average number of random shortcuts per node p. For 

comparison with the previous two extreme cases, we show results for the same composite ratio 9=1/2. The results 
show that for any nonzero value of p, the number of resonances per node will approach a nonzero p > value in the 
large-A^ limit, as opposed to the pure Id ring where it vanishes as 1/N [Fig. (H^a)]. Further, as the number of random 
links per node p increases, p increases monotonically, and the density of resonances initially [0 < p < 0(1)] widens; at 
around p 0(1), the spectrum becomes extended [Fig. [2Kb)]. As we further increase p, the number of resonances per 
node continues to increase monotonically as a function p, quickly "saturating" to its maximum value p = 1 [Fig.[2Ka)], 
while the density of resonances becomes progressively centered about A=0 [Fig. [2Kb,c)], eventually converging to a 
delta-function (if both p—Kxi, N—kx)). Indeed, one can recall for the complete graph, that the average number of 
resonance per node approaches p = 1, but all frequencies are centered about the same value [l^, [l^. Note that in 
both the low shortcut density [0 < p <C 0(1), Fig. [3Ka,b)] and the high shortcut density [p 0(1) (not shown)] 
regimes, for fixed p, the density of resonances becomes independent of the size of the network for large N. Finite-size 
effects are very strong, however, for p ^ 0(1), in particular in the low- and high-frequency regime Fig. [3^0)]. Further 
analysis in the high-connectivity [p ^ 0(1)] regime also reveals that the limit density of resonances has the scaling 
form p(A) = p^/^(/)(p^/^A) [Fig. [21(c), and inset] . This scaling form, valid for all p ^ 0(1), is identical to the one found 
for regular long-range connectivity graphs [ly], including the limit of complete graph (p^N) [29| . 

Next, we provide more details for the low shortcut-density regime (also referred to as the SW regime), < p <C 0(1). 
The scaling of the number of resonances per node p in this regime [Fig. [IJa)], can be extracted from the finite-size 



behavior, typical in Id SW networks [11|, [28|. Since the number of random links per node (density of shortcuts) is 
p, the typical (Euclidean) distance between nodes with shortcuts emanating from them scales as ^ ^ p~^. Thus, for 
N ^ {Np <C 1), there are no random links in almost any realization of the network, and the resonance structure 
will essentially be identical to that of the pure ring. A crossover, governed by the emerging SW structure, can be 
expected when N ^ ^ {Np ^ 1). Thus, in the SW regime, < p <C 0(1), for arbitrary N, the above crossover 
behavior of p(iV, p) can be expressed in terms of N and the scaled variable x = Np with the help of a scaling function 
f{x), such that 

PiN,p) = j^fiNp) , (4) 

where /(a;)~const. for x<^l, while f{x)^x for x^l. Thus, in the large network-size limit (iV— >cxd), for smal/ values of 
p, p ^ p, i.e., the number of resonance per node increases linearly with the average number of random links per node, 
as can be seen in Fig. ^a) (inset). Our numerical results also suggest (although with considerable finite-size effects 
for finite networks) that for any fixed < p <C 1 value, p(A) approaches a system-size-independent limit density and 
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FIG. 3: Density of resonances in Id SW networks for different system sizes, for composite ratio q=l/2, and for (a) p=0.01, (b) 
p=0.10, and (c) p=1.00. System sizes and tlie corresponding symbols in all three panels are the same as in (a). The insets in 
(a) and (b) show the same data on lin-log scales and a Gaussian fit to the largest system size around the center (dashed curve). 



obeys the scaling form p(A) = p{p)p -^/^^'(A/p-^/^) p^/^4'(A/p^/^). Further, the scaling function ^'(s) is reasonably 

2 

well approximated by ~ e^'^" in the vicinity of the center as N^oo, i.e., the density of resonances approaches a 
Gaussian shape [Fig. [3]^a,b)]. 

To model more complicated spatially-embedded random structures [Fig. |D^a)], we considered when both the value 
of the complex link conductivity and the probability to have a link between two nodes can depend on Euclidean 
distance between the two nodes it connects. From elementary length scale considerations, for the distance-dependent 
link conductivity, one has (Jij 1/dij, where dij is the Euclidean distance between nodes i and j (assuming uniform 
"wire" cross-sections) [Eq. ([T]) easily generalizes to this case.] Whereas the probability of having a link (shortcut) be- 
tween node i and j, can also be suppressed e.s^pij ^ ^/dfj (power-law-suppressed SW networks due to "wiring" -cost 
considerations or topological constraints [33,[3H). In Fig.[4jb), we show the resonance spectrum of a two-dimensional 
power-law-suppressed SW network (2d SW) with open boundaries with a=l and p=1.00 (random shortcuts with 
distance-dependent conductivities were added on top of a two-dimensional regular "substrate" [Fig. [4l^a)]), with com- 
posite ratio q—1/2, together with the known results [3, 17 1 of the regular two-dimensional topological structures 
with the same composite disorder. For regular two-dimensional structures, in the large-system size limit, p 0(1) 
and the spectrum is known to be extended 3, 17 1. The addition of distance-dependent shortcuts, however, strongly 
modifies the density of resonances in the vicinity of A=±l (strong peaks for low and high frequencies). Further, the 
structure of the peaks do not approach a limit density in that region, but diverge with system size (with p fixed). 
An analogous plot for an asymmetric link disorder with q=2/3 [Fig. HJ^c)] shows strong (diverging) peaks only in 
the small- frequency regime [also translating to large transient relaxation times or delays in i?C networks Eq. 
Our analyses also indicate that the main qualitative features (articulated peaks for low and/or high frequencies) of 
structures with distance-dependent shortcuts prevail for a range of a, 0<a<ac«2 ± 0.5. 

In summary, we have shown that in random composite networks, by controlling the density of shortcuts p (topological 
randomness) and/or the composite ratio q of the binary links (conductivity disorder), one can effectively shape the 
resonance landscape, or suppress long transient delays in electrical signal propagation. Here, we have highlighted the 
interplay between structural and composite (conductivity) disorder and the collective electrical response in spatially- 
embedded random networks models. The electrical response of more realistic off-lattice random structures, embedded 
in two- and three-dimensions, reflecting relevant wiring cost and topological constraints 1^, 2^, 3^ will be considered 
in future works. A detailed analyses on such structures will help one understand electrical response (resonances and 
signal delays) in complex materials and biological networks. 

This work was supported in part by NSF DMR-0426488. R.H. and S.K.N were also supported in part by the Focus 
Center, NY at RPI. 
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